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ABSTRACT 

We present a systematic examination of the changes in semi-major axis caused by the mutual inter- 
actions of a group of massive bodies orbiting a central star in the presence of eccentricity dissipation. 
For parameters relevant to the oligarchic stage of planet formation, dynamical friction keeps the typ- 
ical eccentricities small and prevents orbit crossing. Interactions at impact parameters greater than 
several Hill radii cause the protoplanets to repel each other; if the impact parameter is instead much 
less than the Hill radius, the protoplanets shift slightly in semi-major axis but remain otherwise un- 
perturbed. If the orbits of two or more protoplanets are separated by less than a Hill radius, they 
are each pushed towards an equilibrium spacing between their neighbors and can exist as a stable 
co-orbital system. In the shear-dominated oligarchic phase of planet formation we show that the feed- 
ing zones contain several oligarchs instead of only one. Growth of the protoplanets in the oligarchic 
phase drives the disk to an equilibrium configuration that depends on the mass ratio of protoplanets 
to planctesimals, S/cr. Early in the oligarchic phase, when E/cr is low, the spacing between rows of 
co-orbital oligarchs are about 5 Hill radii wide, rather than the 10 Hill radii cited in the literature. It 
is likely that at the end of oligarchy the average number of co-orbital oligarchs is greater than unity. 
In the outer solar system this raises the disk mass required to form the ice giants. In the inner solar 
system this lowers the mass of the final oligarchs and requires more giant impacts than previously 
estimated. This result provides additional evidence that Mars is not an untouched leftover from the 
oligarchic phase, but must be composed of several oligarchs assembled through giant impacts. 

Subject headings: planets and satellites: formation — solar system: formation 



1. INTRODUCTION 

The early stages in the formation of planetary 
systems are well described by statistical calculations 
of the evolution of mass distributions and veloc- 
ity dispersions. As larger bodies accumulate from 
the swarm of proto-planetary material, their indi- 
vidual d y namic s begin to dominate their evolution. 
iLissauerl (|1987l ) pointed out that the finite cross- 
section for accretion limits the growth of each proto- 
planet. This is now known as the "oligarchic phase." 
ilKokubo fc Idal Il998f). Numerical dKenvon fc Bromley! 
120061 : iFord fc Chiang) [20071: iLeyison fc Morbidellil l2007t ) 
and analytical (|Goldreich et al.ll2004af ) work has explored 
the transition from oligarchic growth to the chaotic fi- 
nal assembly of the planets. In this work we examine 
the interactions of a moderate number of protoplanets 
in an oligarchic configuration and find that neighboring 
protoplanets stabilize co-orbital systems of two or more 
protoplanets. We present a new picture of oligarchy in 
which each part of the disk is not ruled by one but by 
several protoplanets having almost the same semi-major 
axis. 

Our approach is to systematize the interactions be- 
tween each pair of protoplanets in a disk where a swarm 
of small icy or rocky bodies, the planetesimals, contain 
most of the mass. The planetesimals provide dynamical 
friction that circularizes the orbits of the protoplanets. 
The total mass in planetesimals at this stage is more than 
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that in protoplanets so dynamical friction balances the 
excitations of protoplanets' eccentricities. We character- 
ize the orbital evolution of a protoplanet as a sequence 
of interactions occurring each time it experiences a con- 
junction with another protoplanet. The number density 
of protoplanets is low enough that it is safe to neglect 
interactions between three or more protoplanets. 

To confirm our description of the dynamics and explore 
its application to more realistic proto-planetary situa- 
tions we perform many numerical N-body integrations. 
We use an algorithm optimized for mostly circular orbits 
around a massive central body. As integration variables 
we choose six constants of the motion of an unperturbed 
Keplerian orbit. As the interactions between the other 
bodies in the simulations are typically weak compared 
to the central force, the variables evolve slowly. We em- 
ploy a 4th-order Run ge-Kutta integrati on algorithm with 
adaptive time-steps (|Press et al.l 119921 ) to integrate the 
differential equations. During periods of little interac- 
tion, the slow evolution of our variables permits large 
time-steps. 

During a close encounter, the inter-particle gravita- 
tional attraction becomes comparable to the force from 
the central star. In the limit that the mutual force be- 
tween a pair of particles is much stronger than the central 
force, the motion can be more efficiently described as a 
perturbation of the two-body orbital solution of the bod- 
ies around each other. We choose two new sets of vari- 
ables: one to describe the orbit of the center-of-mass of 
the pair around the central star, and another for relative 
motion of the two interacting objects. These variables 
are evolved under the influence of the remaining parti- 
cles and the central force from the star. 

Dynamical friction, when present in the simulations, 



2 



is included with an analytic term that damps the eccen- 
tricities and inclinations of each body with a specified 
timescale. All of the simulations described in this work 
were performed on Caltech's Division of Geological and 
Planetary Sciences Dell cluster. 

We review of some basic results from the three-body 
problem in section [21 and describe the modifications of 
these results due to eccentricity dissipation. In section 
[31 we generalize the results of the three-body case to 
an arbitrary number of bodies, and show the resulting 
formation and stability of co-orbital sub-systems. Sec- 
tion[4]demonstrates that an oligarchic configuration with 
no initial co-orbital systems can acquire such systems as 
the oligarchs grow. Section O describes our investigation 
into the properties of a co-orbital oligarchy, and section 
[5] places these results in the context of the final stages 
of planet formation. The conclusions are summarized in 
section 

2. THE THREE-BODY PROBLEM 

The circular restricted planar three-body problem 
refers to a system of a zero mass test particle and two 
massive particles on a circular orbit. We call the most 
massive object the star and the other the protoplanet. 
The mass ratio of the protoplanet to the star is /i. Their 
orbit has a semi-major axis a and an orbital frequency f2. 
The test particle follows an initially circular orbit with 
semi-major axis dtp = a(l + x) with x <C 1. Since the 
semi-major axes of the protoplanet and the test particle 
are close, they rarely approach each other. For small x, 
the angular separation between the two bodies changes 
at the rate (3/2)f2a; per unit time. Changes in the ec- 
centricity and semi-major axis of the test particle occur 
only when it reaches conjunction with the protoplanet. 

The natural scale for xa is the Hill radius of the pro- 
toplanet, Rh = (/i/3) 1 / 3 a. For interactions at impact 
parameters larger than about four Hill radii, the ef- 
fects of the protoplanet can be treated as a perturba- 
tion to the Keplerian orbit of the test particle. These 
changes can be calculated analytically. To first order 
in /x, the change in eccentricity is e k = A k [i.x~ 2 , where 
A k = (8/9) [2K (2/3) +#1(2/3)] w 2.24 and K and 
K\ are modified Bessel functions of the second kind 
(|Goldreich fc Tremaind[l978HPetit fc Henonl 1 1986ft . 

The change in semi-major axis of the test particle can 
be calculated from an integral of the motion, the Jacobi 
constant: Cj = E — HH, where E and H and are energy 
and angular momentum per unit mass of the test particle. 
Rewriting Cj in terms of x and e, we find that 

-x 2 — e 2 = const. (1) 
4 

If the encounter increases e, |x| must also increase. The 
change in x resulting from a single interaction on an ini- 
tially circular orbit is 

Ax=(2/3)e 2 k /x=(2/3)A 2 kf i 2 x- 5 . (2) 

The contributions of later conjunctions add to the ec- 
centricity as vectors and do not increase the magnitude 
of the eccentricity by e&. Because of this the semi-major 
axis of the test particle generally does not evolve further 
than the initial change Ax. Two alternatives are if the 
test particle is in resonance with the protoplanet, or if its 



orbit is chaotic. If the test particle is in resonance, the 
eccentricity of the particle varies as it librates. Chaotic 
orbits occur when each excitation is strong enough to 
change the angle of the next conjunctio n substantially; 
in this case e and x evolve stochastically ( Wisdom! 1 19801 : 
iDuncan et aTlll989ft . 

Orbits with x between 2-4 Rh/cl can penetrate the 
Hill sphere and experience large changes in e and a. 
This regime is highly sensitive to initial conditions, so 
we only offer a qualitative description. Particles on 
these orbits tend to receive eccentricities of the order the 
Hill eccentricity, en = Rh/o,, and accordingly change 
their semi- major axes by ~ Rh- We will call this 
the "strong-scattering regime" of separations. A frac- 
tion of these trajectories collide with the protoplanet; 
these orbits are responsible fo r proto-planetary accretio n 
(|Greenzweig fc Lissauer|[l99fl iDones fc Tremain"elll993ft . 

For x < Rh/cl, the small torque from the proto- 
planet is sufficient to cause the particle to pass through 
x = 0. The particle then returns to its original separa- 
tion on the other side of the protoplanet 's orbit. These 
are the famous horseshoe orbits that are related to the 
1:1 mean-motion resonance. The change in eccentricity 
from an initially circular orbit that exp eriences this in- 
teraction can be calculated analytically ([Petit fc Henonl 
[1986^ : e k = 2 2 / 3 3~ 3/2 5r(2/3)^ /3 exp(-(87r/9)/^- 3 ), 
where T(2/3) is the usual Gamma function. Since this in- 
teraction is very slow compared to the orbital period, the 
eccentricity change is exponentially small as the separa- 
tion goes to zero. As in the case of the distant encounters, 
the conservation of the Jacobi constant requires that x in- 
creases as the eccentricity increases (equation [T]) . Then, 

u 2/3 

Ax = 2.83- — exp(-5.58/x:r~ 3 ). (3) 

x 

To apply these results to proto-planetary disks, we 
must allow the test particle to have mass. We now refer 
to both of the bodies as protoplanets, each having mass 
ratios with the central object of fj,\ and /i2- The change 
in their total separation after one conjunction is given by 
equations and [3] with fi — fi\ + \i2- 

Figure [T] plots the change in a after one conjunction 
of two equal mass protoplanets as measured from nu- 
merical integrations. All three types of interactions de- 
scribed above are visible in the appropriate regime of 
x. Each point corresponds to a single integration of two 
bodies on initially circular orbits separated by x. For 
the horseshoe- type interactions, each protoplanet moves 
a distance almost equal to x; we only plot the change 
in separation: Aajj.s. = |Aa| — \x\a. The regimes of 
the three types of interactions are marked in the figure. 
The dashed line in the low x regime plots the analytic 
expression calculated from equation [3) The separations 
that are the most strongly scattered lie between 2 — 4Rh, 
surrounding the impact parameters for which collisions 
occur. For larger separations the numerical calculation 
approaches the limiting expression of equation [2j which 
is plotted as another dashed line. 

The sea of planetesimals modifies the dynamics of 
the protoplanets. If the planetesimals have radii less 
than ~ 1 km, their own collisions balance the excita- 
tions caused by the protoplanets. At the same time, 
the planetesimals provide dynamical friction that damps 



3 



the rate: 



10' 








10 u 


-Horseshoe 


■/ + + + + 




i 






10- 1 


■ / 


+ 




10" 2 


+ / 
+ ■ 


+ 

+ 




x 10- 3 


+ ' 


Strong 




+: 


Rpattprinn 




^ ID' 4 


+:' 
+; 






+; 
+,' 






10- 5 


+■' 
+: 






io- 6 


+ 
+: 

+;' 




Distant - 
Encounters 


10- 7 














IO' 8 









5 

x/R u 



10 



20 



Fig. 1. — The change in semi-major axis after a conjunction 
of two bodies on initially circular orbits whose masses are smaller 
than that of the star by the ratio fj, = 3 X 10 — 9 , plotted as a 
function of the initial separation. The points are calculated with 
numerical integrations, while the dashed lines show the analytic 
results, equations [2] and \3\ At the smallest impact parameters 
the bodies switch orbits; in this case we have measured the change 
relative to the initial semi-major axis of the other protoplanet. The 
horizontal lines separate the regions of x that are referred to in the 
text. 

the eccentricities of the protoplanets. When the typi- 
cal eccentricities of the protoplanets and the planetesi- 
mals are lower than the Hill eccentricity of the proto- 
planets, this configuration is said to be shear-dominated: 
the relative velocity between objects is set by the differ- 
ence in orbital frequency of nearby orbits. In the shear- 
domina ted eccentricity regime, the rate of dynamical fric- 
tion is (jGoldreich et al.ll2004bh : 



1 de 
e dt 



an 2 1 

w — t^ol = — 
pR 



(4) 



where R and p are the radius and density of a proto- 
planet, a is the surface mass density in planetesimals, a 
is the ratio R/ Rh, and Cd is a dimensionless coefficient of 
order unity. Recent studies have found values for Cd be- 
tween 1.2 and 6.2 (Ohtsuki et. al. 2002; Schlichting and 
Sari, in preparation). For this work, we use a value of 1.2. 
For parameters characteristic of the last stages of planet 
formation, Td 3> 27r/f2. The interactions of the proto- 
planets during an encounter are unaffected by dynamical 
friction and produce the change in e and a as described 
above. In between protoplanet conjunctions, the dynam- 
ical friction circularizes the orbits of the protoplanets. 
The next encounter that increases e further increases x 
to conserve the Jacobi constant. The balance between ex- 
citations and dynamical friction keeps the eccentricities 
of the protoplanets bounded and small, but their sepa- 
ration increases after each encounter. This mechanism 
for orbital repulsion has been previously identified by 
iKokubo fc Ida! (|1995f ) . who provide a timescale for this 
process. We alternatively derive the timescale by treat- 
ing the repulsion as a type of migration in semi-major 
axis. The magnitude of the rate depends on the strength 
of the damping; it is maximal if all the eccentricity is 
damped before the next encounter, or -C 4n/(3flx). 
In this case, a protoplanet with a mass ratio p\ and semi- 
major axis a\ interacting with a protoplanet with a mass 
ratio p,2 in the regime of distant encounters is repelled at 



1 da\ Al , , 

-77- = tt^2 (pi +M2 )x 4 fl 

ai at 2tt 



(5) 



For protoplanets in the horseshoe regime, the repulsion of 
each interaction is given by equation^ These encounters 
increase the separation at an exponentially slower rate of: 



dai 
'~dT 



0.67/i 2 (/xi + M 2 )- 2/3 exp(-5.58(/^i + p 2 )x~ 3 )n. 

(6) 

If instead Td » Air/(3Q,x), the eccentricity of the pro- 
toplanet is not completely damped away before the next 
conjunction restores the protoplanet to e ~ e&. The 
rate at which the separation increases is then related 
to the rate of dynamical friction, a oc e-^efx. Qualita- 
tively, this rate is slower than those of equations [5] and [6] 
by (jdSlx)^ 1 . We focus on the maximally damped case 
where Td <C 4ir/(3Qx). 

3. THE DAMPED N-BODY PROBLEM 

Having characterized the interactions between pairs of 
protoplanets, we next examine a disk of protoplanets 
with surface mass density S. Each pair of protoplan- 
ets interacts according to their separations as described 
in section [21 If the typical spacing is of order Rh, the 
closest encounters between protoplanets causes changes 
in semi-major axes of about Rh and eccentricity exci- 
tations to en- The strong scatterings may also cause 
the two protoplanets to collide. If the planetesimals 
are shear-dominated and their mass is greater than the 
mass in protoplanets, the eccentricities of the protoplan- 
ets are held significantly below en by dynamical friction 
(jGoldreich et al.ll2004bl ). and the distribu tion of their ec- 
centricities can be calcul ated analytically (|Collins &: Saril 
120061 : [Collins et al.l2006l ). If the scatterings and collisions 
rearrange the disk such that there are no protoplanets 
with separations of about 2 — ARh, the evolution is subse- 
quently given by only the gen tle pushing of distant inter- 
actions ( Kokubo &: Idal fl995). However, there is another 
channel besides collisions through which the protoplan- 
ets may achieve stability: achieving a semi-major axis 
very near that of another protoplanet. 

A large spacing between two protoplanets ensures they 
will not strongly-scatter each other. However, a very 
small difference in semi-major axis can also provide this 
safety (see figure Q] and equation [S]). Protoplanets sep- 
arated by less than 2Rh provide torques on each other 
during an encounter that switch the order of their semi- 
major axis and reverse their relative angular motion be- 
fore they can get very close. Their mutual interactions 
are also very rare, since their relative orbital frequency 
is proportional to their separation. Protoplanets close 
to co-rotation are almost invisible to each other, how- 
ever these protoplanets experience the same a /a from 
the farther protoplanets as given by equation [5[ We call 
the group of the protoplanets with almost the same semi- 
major axis a "co-orbital group" and use the label N to 
refer to the number of protoplanets it contains. The 
protoplanets within a single group can have any mass, 
although for simplicity in the following discussion we as- 
sume equal masses of each. 
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Different co-orbital groups repel each other at the rate 
of equation El For equally spaced rows of the same num- 
ber of equal mass protoplanets, the migration caused by 
interior groups in the disk exactly cancels the migration 
caused by the exterior groups. We say that the proto- 
planets in this configuration are separated by their "equi- 
librium spacing." We define a quantity, y, to designate 
the distance between a single protoplanet and the posi- 
tion where it would be in equilibrium with the interior 
and exterior groups. The near cancellation of the exte- 
rior and interior repulsions decreases y, pushing displaced 
protoplanets towards their equilibrium spacing. The mi- 
gration rate of a single protoplanet near the equilibrium 
spacing of its group be calculated by expanding equation 
[5] to first order in y and taking the difference between 
interior and exterior contributions: 

-JLM-^aN-^-MlSlN — e H Q, (7) 
y at y f— ^ aix a \ Rh J 

where we assume that the other co-orbital groups in the 
disk are regularly spaced by Act — x a and contain N 
protoplanets of a single mass ratio. Each term in the 
summation represents a pair of neighboring groups for 
which a is evaluated at the unitless separation ix. Since 
the repulsion rate is a sharp function of the separation, 
the nearest neighbors dominate. The coefficient in equa- 
tion [7] takes a value of 121 when only the closest neigh- 
bors are included (i = 1 only). Including an infinite 
number of neighbors increases the coefficent by a factor 
of 1 + 2~ 5 + 3~ 5 + only about 8 percent. 

The dynamics above describe an oligarchic proto- 
planetary disk as a collection of co-orbital groups each 
separated by several Hill radii. It is necessary though to 
constrain such parameters as the typical spacing between 
stable orbits and the relative population of co-orbital sys- 
tems. To determine these quantities we perform full nu- 
merical integrations. Given a set of initial conditions in 
the strong-scattering regime, what is the configuration of 
the protoplanets when they reach a stable state? 

We have simulated an annulus containing twenty pro- 
toplanets, each with a mass ratio of /j, = 1.5 x 1CP 9 to 
the central star. The protoplanets start on circular or- 
bits spaced uniformly in semi-major axis. We dissipate 
the eccentricities of the protoplanets on a timescale of 
80 orbits; for parameters in the terrestrial region of the 
Solar System and using d = 1.2, this corresponds to a 
planetesimal mass surface density of about 8 g cm~ 2 . We 
allow the protoplanets to collide with each other setting 
oT l = 227; this corresponds to a density of 5 g cm -3 . 

We examine two initial compact separations: 1.0 Rh 
(set A) and 2.5 Rh (set B). For each initial separation 
we run 1000 simulations starting from different randomly 
chosen initial phases. After 6 x 10 3 orbital periods the 
orbits of the protoplanets have stabilized and we stop 
the simulations. To determine the configuration of the 
protoplanets, we write an ordered list of the semi-major 
axis of the protoplanets in each simulation. We then 
measure the separation between each adjacent pair of 
protoplanets (defined as a positive quantity) . If the semi- 
major axes of two or more protoplanets are within 2 Rh, 
we assume they are part of the same co-orbital group. 
The average semi-major axis is calculated for each group. 
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Fig. 2. — Histogram of the intra-group and inter-group separa- 
tions between protoplanets in two sets of numerical simulations. 
Each simulation integrates 20 protoplanets with mass ratios of 
3 X 10~ 9 compared to the central mass. They begin on circular 
orbits with uniform separations in semi-major axis; each set of sim- 
ulations consists of 1000 integrations with random initial phases. 
The eccentricities of the protoplanets are damped with a timescale 
of 80 orbits. The smooth line (red) represents the simulations of 
set A, with an initial spacing of 1.0 Rh, an <i the stepped line (blue) 
shows simulations of set B, which have an initial spacing of 2.5 Rh- 

The distance of each member of a group from the aver- 
age semi-major axis we call the "intra-group separation." 
These values can be either positive or negative and, for 
the co-orbital scenarios we are expecting, are typically 
smaller than IRh- 

When one protoplanet is more than 2 Rh from the 
next protoplanet, we assume that the next protoplanet 
is either alone or belongs to the next co-orbital group. 
The spacing between the average semi-major axis of one 
group and the semi-major axis of the next protoplanet or 
co-orbital group we call the "inter-group spacing." These 
separations are by definition positive. 

Finally we create a histogram of both the intra-group 
separations and the inter-group separations of all the 
simulations in the set. For reference, the initial configu- 
ration of the simulations of set B contains no co-orbital 
groups. The resulting histogram would depict no intra- 
group separations, and have only one non-zero bin rep- 
resenting the inter-group separations of x = 2.5Rh- 

Figure [2] shows the histograms of the final spacings of 
the two sets of simulations. The spacings in set A are 
shown in the smooth line (red), and those of set B are 
shown in the stepped line (blue). The initial closely- 
spaced configurations did not survive. The distributions 
plotted in figure [2] reveal that none of the spacings be- 
tween neighboring protoplanets are in the strong scatter- 
ing regime, since it is unstable. This validates the arbi- 
trary choice of 2 Rh as the boundary in the construction 
of figure 2; any choice between 1 and 3 Rh would not 
affect the results. 

The size of the peak of intra-group spacings shows 
that most of the protoplanets in the disk are co-orbital 
with at least one other body. The shape shows that 
the spread of semi-major axis of each co-orbital group 
is small. This is consistent with equation [7j since end- 
point of these simulations is late enough to allow signif- 
icant co-orbital shrinking. The second peak in figure 
represents the inter-group separation. The median inter- 
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group separation in the two sets are A.8Rh and 4.4i?#. 
This is much less than the IQRh usually assumed for the 
spaci ng between protoplan e ts in o ligarchic planet forma- 
tion dKokubo fc Idal fl998l [200l iThommes et~aT1 [200l 
IWeidenschilringll2005h . 

Figure [2] motivates a description of the final configu- 
ration of each simulation as containing a certain number 
of co-orbital groups that are separated from each other 
by 4 — 5Rfj. Each of these co-orbital groups is further 
described by its occupancy number N. For the simula- 
tions of set A, the average occupancy (N) = 2.8, and for 
set B, (N) — 1.8. Since the simulated annulus is small, 
the co-orbital groups that form near the edge are under- 
populated compared to the rest of the disk. For the half 
of the co-orbital groups with semi-major axes closest to 
the center of the annulus, (N) is higher: (N) — 3.5 for 
set A and (N) = 2.0 for set B. 

4. OLIGARCHIC PLANET FORMATION 

The simulations of section[3]demonstrate the transition 
from a disordered swarm of protoplanets to an orderly 
configuration of co-orbital rows each containing several 
protoplanets. The slow accretion of planetesimals onto 
the protoplanets causes an initially stable configuration 
to become unstable. The protoplanets stabilize by reach- 
ing a new configuration with a different average number 
of co-orbital bodies. To demonstrate this process we sim- 
ulate a disk of protoplanets and allow accretion of the 
planetesimals. 

We use initial conditions similar to the current picture 
of a disk with no co-orbital protoplanets, placing twenty 
protoplanets with mass ratios /j = 3 x 10™ 9 on circu- 
lar orbits spaced by 5Rh ■ This spacing is the maximum 
impact para meter at which a proto planet can accrete a 
planetesimal (jGreenberg et al.lll991f ) and a typical stable 
spacing between oligarchic zones (figure [5]). For the ter- 
restrial region around a solar-mass star, this mass ratio 
corresponds to protoplanets of mass 6 x 10 24 g, far below 
the final expected protoplanet mass (see section [6]). Our 
initial configuration has no co-orbital systems. We in- 
clude a mass growth term in the integration to represent 
the accretion of planetesimals onto the protoplanets in 
the regime where the eccentricity of the planetesimals e p 
obeys a 1 / 2 en < e p < ch (jDones k, Tremainelll993l) : 



1 dM _ 2 4 aQ 1 &H 
M dt pR a e p 



(8) 



Protoplanct-protoplanet collisions are allowed. For sim- 
plicity we assume the planetesimal disk does not evolve 
in response to the protoplanets. Eccentricity damping of 
the protoplanets from dynamical friction of the planetes- 
imals is included. The damping timescale, 80 orbits, and 
growth timescale, 4800 orbits, correspond to a planetesi- 
mal surface density of 10 g cm" 2 and a typical planetes- 
imal eccentricity of e p = 5 x 10 ~ 4 . We have again used 
the value Cd = 1.2. These parameters imply a planetes- 
imal radius of ~ 100 m, assuming that the planetesimal 
stirring by the protoplanets is balanced by physical colli- 
sions. Each protoplanet has a density of 5 g cm -3 . The 
annulus of bodies is centered at 1 AU. We simulate 1000 
systems, each beginning with different randomly chosen 
orbital phases. Figure [3] shows the evolution of the semi- 
major axis of the protoplanets in one of the simulations 
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Fig. 3. — Semi-major axes of the protoplanets vs time in a 
simulation of oligarchic growth around a solar mass star. The 
initial mass of each protoplanet is 6 X 10 24 g and each is spaced 
5 Rh from its nearest neighbor. The planetesimals have a surface 
density of 10 g cm -2 and an eccentricity e p = 5 X 10 — 4 . These 
parameters correspond to a damping timescale of 80 years and a 
growth timescale of 4800 years. The sharp vertical lines indicate 
a collision between two bodies; the resulting protoplanet has the 
sum of the masses and a velocity chosen to conserve the linear 
momentum of the parent bodies. 

as a function of time; other simulations behave similarly. 

If there were no accretion, the protoplanets would pre- 
serve their original spacing indefinitely, aside from a slow 
spreading at the edges of the annulus. However, the 
spacing in units of Hill radii decreases as the protoplan- 
ets grow. Eventually their interactions become strong 
enough to cause collisions and large scatterings. This 
epoch of reconfiguration occurs after a time of approxi- 
mately 4000 orbits in the simulation plotted in figure [3l 
At this point the mass of protoplanets has increased by 
roughly a factor of 2.3, meaning the spacing in units of 
Hill radii has decreased by a factor of 1.3. We would 
expect the chaotic reconfiguration to restore the typi- 
cal spacing to about 5Rh by reducing the number of 
oligarchic zones. The figure, in fact, shows 13 zones 
after the first reconfiguration, compared to 20 before. 
Three protoplanets have collided, and four have formed 
co-orbital groups oi N — 2. The co-orbital pairs are vis- 
ibly tightened over the timescale predicted by equation 
[Jj which for the parameters of this simulation is about 
At ks 3 x 10 3 years. The configuration is then stable un- 
til the growth of the bodies again lowers their separation 
into the strong-scattering regime at a time of 1.1 x 10 4 
years. 

The other realizations of this simulation show simi- 
lar results. We find an average co-orbital population of 
(N) = 1.2 in the middle of the annulus after the first 
reconfiguration. This value is lower than those found 
in section [3] because the protoplanets begin to strongly- 
scatter each other when they are just closer than the 
stable spacing. Only a few protoplanets can collide or 
join a co-orbital group before the disk becomes stable 
again. As described in the paradigm of Ko kubo fc Idal 
(1995), a realistic proto-planetary disk in the oligarchic 
phases experiences many such epochs of instability as the 
oligarchs grow to their final sizes. 

5. THE EQUILIBRIUM CO-ORBITAL NUMBER 



G 



As the protoplanets evolve, they experience many 
epochs of reconfiguration that change the typical co- 
orbital number. The examples given in previous sections 
of this work show the result of a single reconfiguration. 
Our choices of initial conditions with the initial co-orbital 
number (N)i = 1 have resulted in a higher final co-orbital 
number (N)f. If instead, (N)i is very high, the final co- 
orbital number must decrease. As the disk evolves, (N) 
is driven to an equilibrium value where each reconfigu- 
ration leaves (N) unchanged. This value, (N) cq , is the 
number that is physically relevant to the proto-planetary 
disk. 

We use a series of simulations to determine (N) eq at a 
fixed value of E and a. Each individual simulation con- 
tains forty co-orbital groups separated by 4 Rjj- This 
spacing ensures each simulation experiences a chaotic re- 
configuration. The number of oligarchs in each group 
is chosen randomly to achieve the desired (N)i. All 
oligarchs begin with e = en and i = ijj to avoid 
the maximal collision rate that occurs if e < a 1/ ' 2 ejj 
(jGoldreich et a l. 2004b). The initial orbital phase, longi- 
tude of periapse, and line of nodes are chosen randomly. 
We set a lower limit to the allowed inclination to prevent 
it from being damped to unreasonably small values. The 
results of the simulations are insensitive to the value of 
this limit if it is smaller than in; we choose 10 -3 in- 

We include an additional force in the simulations to 
prevent the initial annulus from increasing in width. This 
extra force pushes the semi-major axis of a protoplanet 
back into the annulus at a specified timescale. We choose 
this timescale to be longer than the typical time between 
encounters, (fix) -1 , so that multiple protoplanets are not 
pushed to the boundary of the annulus without having 
the chance to encounter a protoplanet a few Hill radii 
away. Collisions between protoplanets are allowed, but 
the protoplanets are not allowed to accrete the planetes- 
imals. Each simulation is stopped when there has not 
been a close encounter for 1.6 x 10 4 orbits. Inspection 
of the simulation results reveals that this stopping crite- 
ria is sufficient for the disk to have reached an oligarchic 
state. We measure the final semi-major axes of the pro- 
toplanets to determine N for each co-orbital group. For 
each set of parameters (E, a, and (N)i) we perform 100 
simulations. 

The numerical values we have chosen for these simu- 
lations reflect planet formation in the terrestrial region. 
We center the annulus of the simulations at 1 AU. We 
adopt the minimum mass solar nebula f or total mass o f 
solids in the annulus, E + cr = 10 g cm" " 2 (lHavashilll98lh . 
and keep this value fixed throughout all the simulations. 
Figure [4] plots the results of simulations for E/cr = 1/10. 
The points connected by the solid line show the average 
(N) f of each set of simulations, while the dashed lines 
show the average value plus and minus one standard de- 
viation of those measurements. For reference, we plot 
another solid line corresponding to (N)i — (N)f. The 
points at low (N)i show a similarity to the results of 
the simulations of sections [3] and [4j stability is reached 
by increasing the number of oligarchs in each co-orbital 
groups. Once (N}i is too high, the chaotic reconfigu- 
ration results in an oligarchy with lower {N). Figure 
H] depicts a feedback cycle that drives (N) towards an 
equilibrium value that remains unchanged by a reconfig- 
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Fig. 4. — The final (N) of simulations against the initial (N) 
for E = 0.9 g cm -2 and <r = 9.1 g cm -2 . For each value of (N)i 
the mass of each protoplanet is adjusted to keep E constant. The 
dashed lines denote the average value plus and minus one standard 
deviation of the measurements. The solid line illustrates where 
{N)i = Wf- 
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Fig. 5. — The final average mass ratio, (fi), of the protoplanets 
plotted against the final ( N ) for ratio of surface densities of E/cr = 
1/10. Each symbol corresponds to a value of (N)i. The solid lines 
plot lines of constant E for values of (x) one standard deviation 
away from the best fit curve of constant E to the simulations with 
(N)i = 2.5. 

uration. For E/cr = 1/10, we find (A) oq w 2.5. The 
intersection of the dotted lines with (N)i = (N)f yields 
the one standard deviation range of (iV) eq , 2 — 3.2. 

The cause of the wide distribution of each (N) / is ev- 
ident from figure [5] In this figure we plot the values of 
(N) f against the average mass of each protoplanet in 
the same simulations of E/cr = 1/10. All of the points 
lie near a single line of (N) t oc (/i}~ 2 / 3 . This relation is 
derived from the definition E = Nm p / (27rAaa). We find 
the relation 

where we have defined xh to be dimensionless and equal 
to Aa/Rff. While the points in figure \5\ generally follow 
the function given by equation[51 there is significant scat- 
ter. We interpret this variation as a distribution of the 
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Fig. 6. — The equilibrium average co-orbital number (JV) eq 
plotted against the surface mass density ratio of protoplanets to 
planetesimals, E/cr. The error bars represent the standard devia- 
tion of (A r ) C q as defined in the text. The solid and dashed points 
correspond to simulations at 1 AU and 25 AU respectively. The 
dashed points are offset by 5 % in E/cr to distinguish them from 
the solid points. 



average spacing between rows, (xjj)f- For the (iV), = 2.5 
simulations, we measure an average {xu)f — 5-4, with a 
standard deviation of 0.2. The solid lines in figure [5] cor- 
respond to the lower and upper bounds of (xh) / given 
by one standard deviation from the mean. This reaffirms 
our earlier conclusion that the spacing between rows is 
an order unity number of Hill radii of an average size 
body. 

The ratio of E/cr increases as the oligarchs accrete the 
planetesimals. To demonstrate the evolution of (iV) eq 
and (xH)eq, we performed more simulations with values 
of E/cr in the range 0.001-2. At each value we examine a 
range of (AT), to determine (N) eq . We plot the resulting 
values in figure [6] The error bars on the points show 
where one standard deviation above and below (N) f is 
equal to (N)i. As the disk evolves and E/cr approaches 
unity, (N) cq decreases. For high values of E/cr, the equi- 
librium co-orbital number asymptotes towards its mini- 
mum value by definition, 1. 

For the simulations with (A% q , we also measure the 
average spacing between co-orbital groups directly. The 
average spacing in units of the Hill radii of the average 
mass protoplanet, (xff} eq is plotted against (N) eq in fig- 
ure[7l Early in the disk, when E/cr is very small, {xu) e q 
is approximately constant at a value of 5.5. The average 
spacing grows however as E/cr approaches unity. 

Figure [5] shows that all oligarchies of a fixed E exhibit 
similar average spacings (xh)- The points from simu- 
lations of different (N)i confirm that a broad range of 
(N) and (/z) can be achieved, with the relation between 
(N) and (/i) given by equation [9] By finding the equilib- 
rium {N} reached by the disk after many configurations, 
we also fix the average mass of the protoplanet, denoted 
(n) C q- We plot (^) eq //XEarth as a function of E/cr at a = 1 
AU in figure0 where fJ^Em-th is the mass ratio of the Earth 
to the Sun. The error bars show the standard deviation 
of (fi) for the simulations with (N)i — (N) cq . 

For comparison, we also plot (/i) as given by equation[9] 
for a constant (N)i — 1 and (xh) = 5. These parameters 
reflect the typical oligarchic picture with no co-orbital 
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Fig. 7. — The equilibrium average spacing between co-orbital 
groups, (x[i) cq in for simulations with {N)i = (N) cq plotted 
against the surface mass density ratio E/cr. The error bars re- 
flect the standard deviation of the measurements of (xjj) of each 
simulation. 




1e-07 



1e-08 



Fig. 8. — The average mass of the protoplanets in an equilibrium 
oligarchy as a function of the surface mass density ratio E/cr at 
a = 1 AU. The error bars are the standard deviation in average 
mass of the simulations for E/cr and (N)i = (N) eq . The solid 
line plots the average protoplanet mass given by an {N) = 1 and 
{xh) = 5 oligarchy commonly assumed in the literature, described 
by equation [9] 

oligarchs and a fixed s pacing in Hill units (lLissauerj|1987l ; 
IKokubo fc Idal [19951: iGoldreich etldl I2004an . At low 
E/cr, the solid line over-estimates the protoplanet mass 
by over an order of magnitude. This is a result of large 
(N) eq , which allows the disk mass to be distributed into 
several smaller bodies instead of a single protoplanet in 
each oligarchic zone. For E/cr greater than about 0.5, the 
lines cross, and the simple picture is an underestimate of 
(^} cq . Although (iV) eq is close to one for these disks, 
(xH)eq grows, increasing the relative amount of the total 
disk mass that has been accreted into each protoplanet. 

We performed the same calculations for several sets of 
simulations with the annulus of protoplanets centered at 
25 AU. The values of {N) oq we find for these simulations 
are plotted as the dashed line in figure [6] For E/cr < 0.1, 
the co-orbital groups tend to contain more oligarchs at 
25 AU than at 1 AU, but the spacing between rows is 
still (xH)eq ~ 5.5. For larger E/cr, the distance of the 



protoplanets from the star matters less. 

6. ISOLATION 

Oligarchic growth ends when the protoplanets have ac- 
creted most of the mass in their feeding zones and the 
remaining planctcsimals can no longer damp the eccen- 
tricities of the protoplanets. The eccentricities of the 
protoplanets then grow unchecked; this is known as the 
"isolation" phase. The mass of a protoplanet at this 
point is referred to as the "isolation mass," and can be 
found from equation [5] 
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E/cr \ Aidisk (xh) 
E/a + lJ M star (N) 
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The literature typically assumes that at isolation all of 
the mass is in protoplanets. This is equivalent to the 
limit of E/cr > 1. 

The results of section [5] show that oligarchy at a fixed 
semi-major axis is uniquely described by E/cr. For the 
terrestrial region then, M- mo is given by the parameters 
we calculate in section [SJ and is plotted as a function of 
E/cr in figure [H 

The exact ratio of mass in protoplanets to that in 
planetesimals that allows the onset of this instability in 
the terrestrial region is not known; simulations suggest 
that in the outer sola r system this fraction E/cr ~ 10 
(|Ford fe; Ch iang 2007). It is not straightforward to de- 
termine the value of E/cr for which isolation occurs. In 
many of our simulations, the eccentricities of the proto- 
planets rise above en, yet an equilibrium is eventually 
reached. We postpone a detailed investigation of the dy- 
namics of the isolation phase for a later work. For any 
value of E/cr at isolation however, the properties of the 
oligarchy at this stage can be read from figures I6|7l and 

m 

The fate of the protoplanets after isolation depends 
on their distance from the star. In the outer parts of 
the solar system, the nascent ice giants are excited to 
high ecce ntricities and may be e j ected from the system 
entirely {G oldreic h et all l2004at iFord fc Chiarig] 120071 : 
iLevison "fc Morbidelli 2007). Their lower rate of col- 
lisions also likely increases their equilibrium co-orbital 
number for a fixed E/cr relative to this work performed in 
the terrestrial region. In contrast to giant impacts, ejec- 
tions do not change the mass of individual protoplanets, 
so they must reach their full planetary mass as oligarchs. 
For an (N) / 1 at isolation, the mass of the disk needs 
to be augmented proportionally to (N) so that (n) eq at 
isolation is equal to the mass of an ice giant. 

The terrestrial planets tend to collide before they 
can be ejected, as the escape velocity from their sur- 
faces is smaller than the velocity needed to unbind 
them from solar orbits ([Chambersll2001t iGoldreich et al.l 
l2004allKenW & Bromlevll2006lh This process conserves 
the total mass of protoplanets so Md; s k is given by the 
Minimum Mass Solar Nebula. Accounting for (N) ^ 1 
in this case reduces the mass of each body at isolation 
proportionally to (N) 3 ^ 2 . This in turn increases the num- 
ber of giant impacts necessary to assemble the terrestrial 
planets. 

7. CONCLUSIONS AND DISCUSSION 



We have analyzed the interactions of a disk of proto- 
planets experiencing dynamical friction. Conjunctions of 
a pair of protoplanets separated by more than 3 Rr in- 
crease the separation of that pair. The repulsions from 
internal protoplanets cancel those from external proto- 
planets at a specific equilibrium semi-major axis. Several 
bodies can inhabit this semi-major axis on horseshoe- 
like orbits. We have shown through numerical simu- 
lations that these co-orbital systems do form and sur- 
vive. We expect the oligarchic phase of planet formation 
to proceed with a substantial population of co-orbital 
protoplanets. We present an empirical relation between 
the ratio of masses in protoplanets and planetesimals, 
E/cr, and the equilibrium average co-orbital number (N) 
and the equilibrium average spacing between co-orbital 
groups (xh)- To form the extra ice giants that popu- 
late the co-orbital groups in the outer solar system, the 
mass of the proto-planetary disk must be enhanced by 
(N) relative to the existing N = 1 picture. To form the 
terrestrial planets requires (iV) 3 / 2 more giant impacts. 
While we have not calculated the critical value of E/cr 
that initiates the isolation phase, we have completely de- 
termined the parameters of a shear-dominated oligarchy 
of protoplanets up to that point. 

In section |3l we have ignored the repulsive distant 
interactions between a protopla net and the planetesi- 
mals t hat cause typ e I migration ([Goldreich fc Tremaind 
119801 IWardl Il986f i. The additional motion in semi- 
major axis is only a mild change to the dynamics. In 
a uniform disk of planetesimals, an oligarchic configu- 
ration of protoplanets migrates inward at an average 
rate specified by the typical mass of the protoplanets. 
Mass variation between the protoplanets of different co- 
orbital groups causes a differential migration relative 
to the migration of the entire configuration. However, 
the repulsion of the neighboring co-orbital groups coun- 
teracts the relative migration by displacing the equi- 
librium position between two groups by an amount ~ 
(cr /T,)(Rh l 'o)Rh ■ Differential migration also acts on 
members of a single co-orbital group, however its effects 
cannot accumulate due to the horseshoe-like co-orbital 
motion. The ratio of the timescale for migration across 
the co-orbital group to the interaction timescale sets a 
minimum safe distance from the equilibrium separation: 
Vs^c/Rh ~ Ai" 1/6 (Misk/A^ ) 1/2 - For typical co-orbital 
group, where y ~ Rh, the migration is never fast enough 
for a protoplanet to escape the group before the next en- 
counter with a co-orbiting protoplanet brings it to the 
other side of the nominal equilibrium semi-major axis. 

It is also possible that the disk of planetesimals is not 
uniform. The accretional growth of a protoplanet may 
lower the surface density of planetesimals at that semi- 
major axis such that the total mass is locally conserved. 
One might naively expect that the deficit of planetesi- 
mals exactly cancels the repulsion caused by the formed 
protoplanet. However, it can be seen from equation [5] 
that the rate of repulsion of a protoplanet from another 
protoplanet of comparable mass is twice that of the same 
mass in planetesimals. The net rates of repulsion of the 
protoplanets in this scenario are reduced by a factor of 
two; the dynamics are otherwise unchanged. 

One important question is that of the boundary con- 
ditions of a planet-forming disk. The initial conditions 
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of the simulations we present only populate a small an- 
nulus around the central star. We artificially confine the 
bodies in this region to force the surface mass density to 
remain constant. The behavior of S over a larger region 
of the disk may not be similar to that of our annulus. The 
presence of gas giants or previously formed planets may 
prevent any wide-scale diffusion of protoplanets across 
the disk. On the other hand, the dynamics in a logarith- 
mic interval of semi-major axis may not be affected by 
the populations internal and exterior to that region. The 
behavior of protoplanets in the oligarchic phase in a full 
size proto-planetary disk is an open question. 

Earlier analytical work has examined the interac- 
tions between oligarch s that share a feeding zone 
(|Goldreich et al.l [20 04b). These authors conclude that 
protoplanets in an oligarchic configuration are always re- 
duced to an (N) = 1 state. However, we have shown 
that for a shear-dominated disk, the collision rate be- 
tween protoplanets is suppressed as the protoplanets are 
pushed towards almost the same semi-major axis. The 
growth rate of the protoplanets of each co-orbital group 
depends on the eccentricity of the planetesimals. For 
e, p < a 1 / 2 eH the growth rate of a protoplanet scales as 
Br 1 . This is called "orderly" growth since all of the 
protoplanets approach the same size. In the intermedi- 
ate shear- dominated regime of a 1 / 2 en < e P < e#, the 
growth rate is independent of R. The protoplanets then 
retain the relative difference in their sizes as they grow. 
For shear-dominated disks, which are the focus of this 
paper, the co-orbital groups are not disrupted by differ- 
ential growth. 

The spacing between co-orbital groups that we observe 
for most E/<7 is smaller than th e 10 Rh that is typically 
assumed (IKokubo fc Idal 1 19981 l2002t iThommes et all 
2003; Weidenschilling 2005) based on the simulations 
bv iKokubo fc Idal (|1998f ). Their simulations are m 
the dispersion-dominated eccentricity regime, where the 
maximum distance at which an oligarch can accrete a 
planetesimal is set by the epicyclic motion of the plan- 
etesimals, ~ ea. This mo tion sets the width of the feed- 
ing zones; the figures of IKokubo fcldal (fl998l) indicate 
that the typical eccentricity of the smaller bodies cor- 
responds to a distance of IQRh- Dispersion-dominated 
disks with different values for protoplanet sizes and plan- 
etesimal eccentricities should undergo oligarchy with a 
different spacing. In shear-dominated disks, we have 
shown that separations of about 5Rh are set by the dis- 
tant encounters with the smallest impact p arameters. 

The simulations of IKokubo fc I da (1998) do not con- 
tain any co-orbital groups of protoplanets; this is ex- 
pected due to the small number of protoplanets that 
form in their an nulus and the fac t that their eccentricities 
are super- Hill. IThommes et~aLl (|2003f) examine a broad 
range of parameters of oligarchic growth, but the number 
of planetesimals are not enough to damp the protoplanet 
eccentricities sufficiently. However, upon inspection of 
their figure 17 we find hints of the formation of co-orbital 
groups. Also, even though a range of separations are vis- 
ible, many adjacent feeding zones are separated by only 
5Rh as we are finding in our simulations. 

Simulations of the oligarchic phase and the isolation 
epoch that follows by iFord fc ChianeJ (|2007h include five 
bodies that are spaced safely by 5Rh- We would not 



expect the formation of co-orbital oligarchs from an ini- 
tial st ate of so few. Interestingly, iLevison fc Morbidellil 
( 2007) use a population of "tracer particles" to calculate 
the effects of planetesimals on their protoplanets and find 
a strong tendency for these objects to cluster both in co- 
orbital resonances with the protoplanets and in narrow 
rings between the protoplanet orbits. This behavior can 
be understood in light of our equation[2]with the dynam- 
ical friction of our simulations replaced by the collisional 
damping of the tracer particles. 

Simulations of moderate numbers of protoplanets with 
eccentricity damping and forced semi-ma j or ax is migra- 
tion were studied bv ICresswell fc Nelson! (|2006l ): indeed 
they observe many examples of the co-orbital systems 
we have described. We offer the following comparison 
between their simulations and this work. Their migra- 
tion serves the same purpose as the growth we included 
in the simulations of section [4] namely to decrease the 
separations between bodies until strong interactions re- 
arrange the system with stable spacings. The co-orbital 
systems in their simulation likely form in the same way 
as we have described: a chance scattering to almost the 
same semi-major axis as another protoplanet. They at- 
tribute the tightening of their orbits to interactions with 
the gas disk that dissipates their eccentricity, however, 
this is unlikely. Although very close in semi-major axis, 
in inertial space the co-orbital protoplanets are separated 
by ~ a for most of their relative orbit. Since the tight- 
ening of each horseshoe occurs over only a few relative 
orbits, it must be attributed to the encounters with the 
other protoplanets, which occur more often than the en- 
counters between the co-orbital pairs. 

Cresswell and Nelson also find that their co-orbital 
pairs settle all the way to their mutual L4 and L5 La- 
grange points; the systems that we describe do not. In 
our simulations a single interaction between neighbors 
moves each protoplanet a distance on the order of the 
width of the largest possible tadpole orbit, Aa/a ~ /i 1 / 2 . 
The objects in the simulations by Cresswell and Nelson 
have much larger mass ratios with the central star and 
larger separations. In their case a single interaction is not 
strong enough to perturb the protoplanets away from the 
tadpole-like orbits around the Lagrange points. We have 
performed several test integrations with parameters sim- 
ilar to those run by Cresswell and Nelson and confirmed 
the formation of tadpole orbits. Finally, their simulations 
model the end of the planet formation and hint at the 
possibility of discovering extrasolar planets in co-orbital 
resonances. In a gas depleted region, we do not expect 
the co-orbital systems that form during oligarchic growth 
to survive the chaos following isolation. 

In the terrestrial region of the solar system, geological 
measurements inform our understanding of the oligarchic 
growth phase. Isotopic abundances of the Martian me- 
teorites, in particular that of the Hafnium (Hf) to Tung- 
sten (W) radioactive system, depend on the timescale for 
a planet to separate internally into a core and mantle . 
Based on these measurements, lHallidav fc Kleinel (2006) 
calculate that Mars differentiated quickly compared to 
the timescale of the Hf-W decay, 9 Myrs. The oligarchic 
picture of equation [9] with (N) = 1 shows that at 1.5 AU 
with (N) = 1, and E ~ a, (/i) w Mmdxs/Mq; accordingly 
these authors infer that Mars was fully assembled by the 
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end of the oligarchic phase and did not participate in 
the giant impacts that assembled Earth and Venus. A 
co-orbital oligarchy, however, lowers the mass of each 
protoplanet at isolation by a factor of (7V) 3//2 . In this 
picture Mars formed through several giant impacts. This 
scenario is consistent with the isotopic data if Mars can 
experience several collisions in 10 Myrs; the collisional 
timescales for (N) > 1 systems merit further investiga- 
tion. 

The rate and direction of the rotation of Mars, how- 
ever, pr ovide further eviden c e for a history of giant im- 
pacts. iDones fc Tremainel (|1993l ) calculate the angu- 
lar momentum provided by the collision-less accretion 
of planetesimals and show that, for any planetesimal 
velocity dispersion, this process is insufficient to pro- 
duce the observed spins. The moderate prograde rota- 
tion of Mars is thus inconsi s tent w ith pure accretionary 
growth. ISchlichting fc Saril ()2006l ) show that the colli- 
sions of planetesimals inside the Hill sphere as they ac- 
crete produces protoplanets that are maximally rotating, 
which is still inconsistent with the current rotation of 
Mars. Giant impacts later re-distribute the spin-angular- 
momentum of the protoplanets but with a prograde bias; 
this then implies that Mars did participate in the gi- 
ant impact phases of the terrestrial region. Again, fur- 
ther studies are necessary to characterize the timescale 
of the collisional period following the isolation phase in 
an (N) > 1 scenario. 

The compositions of the planets offer more clues to 
their formation. As protoplanets are built up from 
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smaller objects in the proto-planetary disk, their com- 
position approaches the average of the ma terial from 
which they accrete. Numerical simulations by Chambers 
(|2001h show that the collisional assembly of protoplanets 
through a {N) = 1 oligarchy mixes material from a wide 
range of semi-major axes. The composition of the planets 
then reflects some average of all available material. The 
three stable isotopes of oxygen are thought to be initially 
heterogeneous across the proto-planetary disk, and offer 
a measurable probe of compositional differences between 
solar system bodies. In the case of the Earth and Mars, 
a small but finite difference in the ratios of these isotopes 
is usually attrib uted to the statistical fluctuations of th e 
mixing process (|Franchi et al.ll2001b IcTzima et al.ll2007l) . 
An (N) > 1 oligarchy requires more collisions; the same 
isotopic variance between Earth and Mars may require a 
larger dispersion in the composition of the smallest proto- 
planetary materials. However, it is necessary to deter- 
mine the extent of spatial mixing in the (N) > 1 picture 
and to understand the cha nges in composition resulting 
from a single giant impact (|Pahlevan fc Stevenson|[2007l) 
before we can estimate the primordial compositional vari- 
ations allowed by this model. 
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